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PARTICLE  SIZE  AND  SURFACE  AREA  ANALYSIS  OF  AMORPHOUS  AND 

CRYSTALLINE  BORON 


ABSTRACT 

Surface  areas  were  measured  and  scanning  electron  microscopy  (SEM)  photographs  were  obtained 
from  seven  boron  particle  samples  used  in  recent  boron  particle  combustion  experiments.  Some  significant 
differences  were  seen  between  the  manufacturer’s  claimed  particle  sizes  and  the  distributions  obtained  from 
the  SEM  photographs.  In  those  samples  with  large  particle  size  distributions,  calculated  surface  areas  differed 
significantly  from  measured  surface  areas  and  it  is  clear  that  a  single  particle  diameter  parameter  is  inadequate 
to  describe  the  corresponding  surface  area  of  the  true  particle  size  distribution.  In  general,  the  particle  size 
distribution  estimates  taken  from  SEM  yield  lower  surface  areas  than  were  measured.  Further  data  on  the 
purity  of  the  samples  are  necessary  to  determine  whether  the  phase  or  purity  of  the  boron  particle  has  an  affect 
on  the  efficiency  of  boron  particle  combustion. 

INTRODUCTION 

Solid  boron  is  a  dense,  energetic  material  with  a  large  heat  of  fusion  and  has  been  under  scrutiny  for 
over  30  years  as  a  potential  solid  or  liquid  propellant  additive.  The  potential  energy  release  from  boron- 
containing  fuels  is  about  36%  greater  in  cal/g  than  the  liquid  hydrocarbon  fuel  JP5  and  86%  greater  than  solid 
fuels  containing  metallic  aluminum.  When  calculated  on  a  volumetric  basis  in  units  of  cal/cm3,  boron  is  3.9 
times  more  energetic  than  JP5.1 

Boron  combustion  occurs  in  two  stages.  Due  to  the  high  melting  point  of  solid  boron  (2400  K)  and 
the  equally  high  vaporization  point  of  B203,  the  first  stage  of  boron  combustion  is  inhibited  by  a  thin  film  of 
liquid  B203  that  forms  on  the  outer  surface  of  the  boron  particle.  This  film  remains  until  temperatures  greater 
than  2400  K  are  reached,  the  film  vaporizes,  and  combustion  proceeds  more  rapidly.  The  second  stage  of 
boron  combustion  begins  when  the  pure  boron  material  freely  combusts  at  the  higher  temperatures.  Several 
groups  have  measured  ignition  of  boron  particles  under  various  conditions.  Macek  et  al.  3  carried  out  work 
in  a  flat  flame  burner  at  medium  pressures  using  laser  ignition.  Li  and  Williams 4,5  examined  small  particles 
in  a  flat  flame  burner  at  low  pressure.  Yuasa  and  Isoda6  obtained  spatially  resolved  spectra  of  flash-ignited 
boron  particles  in  a  stagnation  stream.  Expansion  into  the  high  pressure  regime  has  been  carried  out  by  Kner, 
Foelsche,  and  coworkers.  Krier  et  al. 7  examined  the  effect  of  pressure  on  ignition  and  combustion  in  a  shock 
tube  apparatus.  Foelsche  et  al.8  examined  the  combustion  of  boron  in  hydrogen/oxygen/water  vapor  mixtures 
at  1360  atm  and  3800  K  in  a  combustion  bomb. 

There  is  some  controversy  over  the  details  of  the  first  and  second  stage  and  how  the  B2O3(0  layer  is 
formed.  Based  on  theoretical  calculations,  Glassman9  concluded  that  in  stage  1,  the  boron  oxide  is  formed 
by  the  dissolution  of  pure  boron  outward  through  the  B203(Z)  layer.  However,  King10 12  modeled  the  problem 
by  assuming  02  diffuses  inward  through  the  Bp/O  to  reach  the  B/  B  P3  interface.  Recent  measurements  and 
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modeling  by  Yeh  and  Kuo13  support  the  former  description.  Using  environmental  scanning  electron 
microscopy  (ESEM),  nonconductive  surfaces  were  detected.  With  this  they  observed  the  melting  of  boron 
particles  at  a  temperature  of  950  K,  far  lower  than  the  melting  point  of  boron  and  higher  than  the  melting  point 
of  B203.  They  observed  a  polymeric  vitreous  (BO)n  complex  formed  by  the  reaction  of  dissolved  boron  and 
molten  B2O3(0,  and  report  the  dependency  of  first-stage  combustion  times  on  ambient  gas  temperature, 
ambient  species  concentration,  and  chamber  pressure.  The  pressures  included  in  the  modeling  were  in  the 
range  of  1  to  35  atm.  Foelsche  et  al.8  support  a  model  in  which  the  diffusion  of  02  into  the  boron/B2O3(0  layer 
occurs  during  the  first  stage.  This  is  based  on  the  fact  that  they  saw  a  two  times  increase  in  ignition  delay  with 
a  reduction  of  5%  in  the  02  mole  fraction.  Their  data  were  collected  under  conditions  of  high  temperature 
(1400  to  3400  K)  and  pressure  (peak  of  157  atm). 

During  combustion,  a  variety  of  gas-phase  products  are  produced  including  but  not  limited  to  BO, 
B02,  B202,  BOH,  and  HOBO.14'16  The  lowest  energy,  and  therefore,  the  most  desired  product  is  B203.  The 
reaction  from  solid  boron  to  liquid  B203  is  exothermic  by  300  Kcal/mol  and  is  shown  in  Eq.  (1). 


2  B(s)  +  3/2  02  -B2O3(<0  (1) 


This  is  the  most  efficient  pathway  in  the  combustion  of  boron.  The  more  realistic  and  undesired  scenario  is 
to  form  a  variety  of  species  other  than  B203,  for  example,  HOBO.  The  gas-phase  modeling  of  boron 
combustion  is  based  on  chemical  kinetics  code  packages  developed  at  Sandia  called  CHEMKIN17  and 
SENKEN.18  Fluorine  chemistry  has  been  added  recently  by  Yetter  et  al.14  and  Brown  et  al.15  to  demonstrate 
that  enhancement  of  complete  combustion  can  occur  when  HF  and  other  fluorinated  compounds  are  added 
to  the  pool  of  reactants.  A  few  other  approaches  have  been  investigated  to  try  to  enhance  the  performance  of 
boron  combustion.  One  of  these  is  to  add  a  thin  layer  of  titanium  or  other  metal  to  the  surface  of  the  boron 
particle  to  heat  the  surface  faster,  and  hence,  reach  the  second  stage  of  combustion  faster.19,20 

It  has  been  suggested  that  one  difference  in  boron  particle  combustion  performance  may  be  the  phase 
and  related  purity  of  the  boron  material.  In  general,  it  has  been  difficult  to  generate  amorphous  material  with 
as  high  a  purity  as  crystalline  material.  Whether  there  is  a  phase/purity  related  difference  in  performance 
remains  controversial  because  few  controlled  experiments  have  been  done  to  prove  or  disprove  this 
hypothesis.  Much  of  the  older  literature  and  some  of  the  newer  literature  are  complicated  by  the  fact  that 
boron  samples  are  not  fully  characterized  and  in  cases  where  characterization  is  done,  the  results  are 
conflicting.  Small  impurities  appear  to  strongly  affect  some  of  the  physical  properties  of  boron.21  Literature 
in  the  1960's  and  1970's  is  difficult  to  interpret  but  folklore  suggests  that  contaminants  in  amorphous  samples 
were  the  cause  of  increased  performance.  These  difficult  questions  of  physical  properties  and  purity  cloud 
the  literature  and  remain  unanswered  today. 

In  this  report,  seven  samples  of  amorphous  and  crystalline  boron  that  recently  have  been  used  in 
various  combustion  and  detonation  experiments  are  analyzed  by  SEM  and  surface  area  analysis.  The  boron 
particles  were  examined  by  SEM  to  determine  the  general  physical  features  and  the  particle  size  distribution. 
An  estimated  surface  area  (m2/  gm)  was  calculated  from  the  measured  particle  size  distribution.  The  actual 
surface  area  was  then  experimentally  measured  using  a  Monosorb  instrument  (Quantachrome  Corp.)  and 
analyzed  using  a  single-point  BET  (Brunauer,  Emmett  and  Teller)  analysis.22  The  actual  and  calculated 
surface  areas  are  compared  and  discussed  in  terms  of  recent  particle  combustion  results. 
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Yeh  and  Kuo13  have  theorized  that  if  the  particle  reactions  are  diffusion  controlled,  the  particle 
diameter  has  a  larger  effect  (<*d2)  on  the  combustion  time.  However,  if  kinetically  controlled,  the  combustion 
time  is  linear  with  particle  diameter  («d).  The  two  pathways  are  equal  when  the  Damkohler  number  equals 
one  or  P*dp/75  equals  one,  where  P  is  in  atm  and  d  p  is  in  pm.  Therefore,  the  lower  the  pressure,  the  more 
influence  the  surface  area  will  have.  One  way  to  increase  the  surface  area  is  to  decrease  the  particle  size. 


EXPERIMENT 
Boron  Samples 

Table  1  lists  the  seven  samples  that  were  analyzed  in  this  study.  These  samples  provided  a  variety 
of  surface  areas  from  stocks  that  had  been  used  in  boron  combustion  experiments. 

Surface  Area  Measurements 

The  surface  area  measurements  were  made  using  a  Monosorb  instrument  that  uses  a  single-point  BET 
algorithm.22  A  mixture  of  30%  nitrogen  (by  volume)  with  the  balance  helium  is  used  as  the  adsorbate.  A 
sample  of  about  0.25  to  0.3  grams  was  degassed  at  18°C  for  45  min  prior  to  measuring  the  surface  area  to 
remove  water  and  adsorbed  gases  at  the  surface.  The  sample  was  then  cooled  to  liquid  nitrogen  temperatures 
and  the  adsorbate  adhered  to  the  boron  surface.  Once  at  equilibrium,  the  sample  was  heated  and  the  adsorbate 
was  quickly  desorbed.  The  flow  of  gas  out  of  the  sample  was  measured  with  thermal  conductivity  detectors 
and  the  surface  area  was  calculated  from  the  amount  and  cross  sectional  area  of  the  adsorbate  released.  The 
measurement  was  carried  out  at  a  partial  pressure  of  adsorbate  of  about  0.3.  This  ensured  that  a  monolayer 
had  formed  on  the  surface  and  allowed  the  calculation  of  the  surface  area  based  on  the  cross-sectional  area  of 
the  adsorbate.  The  amount  of  gas  released  was  converted  to  units  of  square  meters.  The  sample  was  then 
weighed  and  the  surface  area  reported  in  units  of  square  meters  per  gram.  The  accuracy  of  the  single-point 
BET  measurement  is  generally  on  the  order  of  5%.  The  reproducibility  is  ~1  to  2%. 


RESULTS 

Scanning  Electron  Microscopy 

Scanning  electron  microscopy  was  carried  out  on  a  Hitachi,  Model  S400  at  15  to  20  kV  at 
magnifications  from  400x  to  30,000x.  Sample  1  was  coated  with  spuddered  gold-palladium  for  improved 
contrast  and  is  shown  in  Fig.  1.  The  particles  are  spherical  in  appearance  and  relatively  uniform  with  dp’s  ~ 
0.08  to  0.1  pm.  Sample  2  is  shown  in  Fig.  2.  It  appears  to  be  identical  in  size  and  structure  to  sample  1. 
Figures.  3(a)  and  3(b)  are  that  of  sample  3,  indicating  a  larger  size  distribution  than  samples  1  and  2.  Figures 
4(a)  and  (b)  represent  the  crystalline  boron  sample  4,  which  is  0.08  to  8.0  pm  in  particle  size  and  appears  as 
a  very  wide  distribution  of  irregularly  shaped  crystals.  Figures  5(a)  and  (b)  represent  a  crystalline  boron 
sample  similar  to  sample  4.  Figures  6(a)  and  6(b)  represent  crystalline  boron  with  dp’s  from  0.08  to  35.0  pm. 
Figures  7(a)  and  7(b)  represent  a  sieved  sample  of  larger  crystalline  boron  particles  with  dp’s  between  20  and 
25  pm.  However,  as  can  be  seen  in  Fig.  7(b),  there  are  numerous  small  particles  which  adhere  to  the  larger 
crystalline  surfaces  which  are  not  removed  by  sieving. 


Table  1 

Particle  Sizes  and  Surface  Areas  of  Amorphous  and  Crystalline  Boron 
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Fig.  1  -Shown  is  a  SEM  photograph  of  sample  1,  Fig.  2-  Shown  is  a  SEM  photograph  of  sample  2, 

amorphous  boron,  from  Zavitsanos.  It  is  coated  with  amorphous  boron  from  Krier,  with  a  magnification  of 

spuddered  gold-palladium  for  improved  contrast.  It  was  30,000x. 

taken  at  15  kV  with  a  magnification  of  30,000x.  A  1  pm 
scale  is  given  by  dots  in  the  lower  right  hand  comer. 
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Fig.  3(a)  Shown  is  a  SEM  photograph  of  sample  3,  Fig.  3(b)-Shown  is  a  SEM  photograph  of  sample  3, 

amorphous  boron  from  Kuo  and  Yeh,  with  a  amorphous  boron  from  Kuo  and  Yeh,  with  a 

magnification  of  4,000x.  magnification  of  1 2,000x. 
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particle  size,  dp(pm) 


particle  size,  dp(pm) 


Fig.  3(c)-The  upper  plot  shows  the  particle  size  distribution  of  sample  3  and  the  lower 
plot  shows  the  surface  area  of  the  total  particle  size  distribution  normalized  by  the 
interval  size.  See  Appendix  A  for  a  tabular  list. 
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Fig.  4(a)-Shown  is  a  SEM  photograph  of  sample  4,  Fig.  4(b)-Shown  is  a  SEM  photograph  of  sample  4, 

crystalline  boron  from  Russell,  with  a  magnification  of  crystalline  boron  from  Russell,  with  a  magnification  of 

3,000x.  15,000x. 


Particle  Size  and  Surface  Area  Analysis  ofAmorphous  and  Crystalline  Boron 


9 


120 


100 

</) 

<1) 

O  80 
2.  60 


d) 

.Q 

E  20 

u 

c 

0 


0  2  4  6  8  10 

particle  size,  dp(pm) 


x"E  3.2e-1 1 
.2  "T  2.8e-11 

O  (0 

t:  ®  2.4e-11 
to  o 

■S-t  2.0e-1 1 
CO  co 

2  1 .6e-1 1 

ro  o 

0)  1 .2e-1 1 

o  <D 

<2  -Q  8.0e-12 

^  E 

«  |  4.0e-12 

0.0 


-2  0  2  4  6  8  10 

particle  size,  dp(pm) 


Fig.  4(c)-The  upper  plot  shows  the  particle  size  distribution  of  sample  4  and  the  lower  plot  shows 
the  surface  area  of  the  total  particle  size  distribution  normalized  by  the  interval  size. 
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Fig.  5(a)-Shown  is  a  SEM  photograph  of  sample  5,  Fig.  5(b)-Shown  is  a  SEM  photograph  of  sample  5, 

crystalline  boron  from  Kuo  and  Yeh,  with  a  magnification  crystalline  boron  from  Kuo  and  Yeh,  with  a  magnification 

of  2,000x.  of  6,000x. 


surface  area/particle 
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particle  size,  dp(pm) 


particle  size,  dp(pm) 


Fig.  5(c)-The  upper  plot  shows  the  particle  size  distribution  of  Sample  5  and  the  lower  plot 
shows  the  surface  area  of  the  total  particle  size  distribution  normalized  by  the  interval  size. 
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Fig.  6(a)-Shown  is  a  SEM  photograph  of  sample  6,  Fig.  6(b)-  Shown  is  a  SEM  photograph  of  sample  6, 

unsieved  crystalline  boron  from  Krier,  with  a  unsieved  crystalline  boron  from  Krier,  with  a 

magnification  of  400x.  magnification  of  6,000x. 


surface  area/particle  x 

number  of  particles,  (m  )  number  of  particles 
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particle  size,  dp(pm) 


Fig.  6(c)-The  upper  plot  shows  the  particle  size  distribution  of  sample  6  and  the  lower  plot  shows 
the  surface  area  of  the  total  particle  size  distribution  normalized  by  the  interval  size. 
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Fig.  7(a)-Shown  is  a  SEM  photograph  of  sample  7,  Fig.  7(b)-Shown  is  a  SEM  photograph  of  sample  7, 

sieved  crystalline  boron  from  Krier,  with  a  magnification  sieved  crystalline  boron  from  Krier,  with  a  magnification 

of  400x.  of  6,000x. 
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Fig.  7(c)-The  upper  plot  shows  the  particle  size  distribution  of  sample  7  and  the  lower 
plot  shows  the  surface  area  of  the  total  particle  size  distribution  normalized  by  the 
interval  size. 
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Estimating  the  Size  Distribution  of  Particles 

Samples  1  and  2  contained  amorphous  particles  that  were  uniform  and  spherical,  therefore,  the  size 
distribution  was  represented  by  a  Delta  function.  The  particle  size  of  highest  occurrence,  dp(max),  was  0.85 
pm.  In  sample  3,  the  particles  were  not  uniform  and  a  distribution  determined  from  SEM  was  required  to 
describe  the  particle  sizes.  For  each  of  the  ciystalline  boron  samples  4  through  6,  the  distribution  of  particles 
was  broad  and  also  required  an  estimate  of  the  distribution  from  SEM  photographs.  Both  the  Delta  function 
and  distribution  methods  were  applied  to  sample  7.  Table  2  lists  the  observed  particle  size  distributions;  the 
10%/90%  particle  size  distribution;  the  most  probable  particle  size,  dp(max);  the  d^  particle  diameter  in  which 
the  surface  area  is  normalized  to  the  number  of  particles;  and  the  d32  particle  diameter,  also  known  as  the 
Sauter  diameter,  in  which  the  total  surface  area  is  normalized  to  the  mass  of  the  particle  ensemble.  Reference 
23  and  Appendix  A  discuss  these  diameters. 

Particle  size  distributions  were  obtained  by  manually  counting  the  particles  in  the  SEM  photographs. 
A  minimum  of  100  particles  was  counted  in  each  sample  to  obtain  the  distributions.  In  sample  3,  larger  sized 
particles  were  counted  from  an  area  of  288  pm2,  where  7.5  pm  equaled  3.1  cm  in  the  photograph.  Smaller 
particles  were  counted  in  an  area  of  30.2  pm2  where  2.5  pm  equaled  3.1  cm.  The  counts  in  these  areas  were 
then  multiplied  by  the  ratios  of  the  areas  to  arrive  at  the  total  particle  distribution  in  the  1022  pm2  area.  The 
particle  diameters  were  measured  in  one  dimension  along  one  axis  of  the  photograph.  Since  the  particles  were 
of  random  shapes,  and  were  randomly  positioned  in  the  photo,  this  simpler  measurement  was  assumed  to  give 
an  average  particle  diameter.  The  particles  were  then  assumed  to  be  spherical,  which  represents  a  minimum 
calculated  surface  area  for  the  particles  observed  since  any  deviation  away  from  a  perfect  sphere  will  lead  to 
a  larger  surface  area. 

In  sample  4,  large  particles  were  counted  in  a  photographed  area  of  29.2  by  35  pm  equaling  1022  pm2, 
where  10  pm  corresponds  30  mm  in  the  photograph.  Mid-sized  particles  were  counted  in  an  area  of  51 1  pm2 
and  small  particles  were  counted  in  an  area  of  41  pm2.  Figure  4(c)  is  a  plot  of  the  particle  distribution  for  this 
sample.  In  sample  5,  large  particles  were  counted  in  an  area  corresponding  to  246  pm2  (l/9th  the  area  of  that 
of  Fig.  5(a))  and  small  particles  were  counted  from  all  of  Fig.  5(b),  also  corresponding  to  246  pm2.  The 
particle  size  distribution  of  sample  6  was  counted  similarly.  An  area  of  10,900  pm2  was  used  to  count  the 
larger  particles  and  an  area  of  254  pm2  was  used  for  the  smaller  particles.  The  total  distribution  of  particles 
for  5450  pm2  is  shown  in  Fig.  6(c).  The  particle  size  distribution  of  sample  7,  sieved  crystalline  boron,  is 
shown  in  Fig  7(c).  The  distribution  of  particles  were  obtained  from  the  SEM  photographs  7a  and  b  in  a  similar 
manner  described  for  Sample  6,  however,  an  area  of  5450  pm2,  not  10,900  pm2,  was  used.  Seventy-five 
microns  corresponds  to  30  mm  in  Fig.  7(a)  and  5  pm  corresponds  to  30  mm  in  Fig.  7(b). 

Estimating  the  Surface  Area  from  the  Size  Distributions  of  the  Particles 

Using  the  particle  size  distributions  obtained  from  SEM  photos,  and  assuming  the  particles  have  a 
spherical  shape,  the  surface  area  of  each  boron  sample  in  units  of  m2/g  was  calculated.  Larger  particles  have 
larger  surface  areas,  however,  smaller  particles  have  larger  surface  areas  per  gram.  A  plot  of  the  surface  area 
per  particle  as  a  function  of  particle  diameter  is  shown  in  the  upper  plot  of  Fig.  8.  This  surface  area  was 
multiplied  by  the  number  of  particles  to  obtain  the  total  surface  area  shown  in  the  lower  plots  of  Figs.  3(c), 
4(c),  5(c),  6(c),  and  7(c).  This  was  related  to  the  number  of  grams  in  the  sample  by  calculating  the  number  of 
grams  of  boron  material  present  in  the  photographs  from  the  particle  diameter,  particle  distribution,  and  the 
density  of  boron  (2.34  g/cm3).  To  calculate  the  surface  area  of  uniformly  sized  particles,  the  following 


surface  area/particle,  (m 
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Fig.  8-The  upper  plot  shows  the  surface  area  of  a  simple  spherical  particle  as  a  function 
of  particle  size  and  the  lower  plot  shows  the  surface  area  per  gram  of  a  simple  sphere  as 
a  function  of  particle  size. 
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equation  derived  from  the  surface  area,  volume  of  a  sphere,  and  density  was  used 

2.56/dp(pm)  =  surface  area  (m2/g).  (2) 

A  general  plot  of  surface  area  per  gram  of  boron  as  a  function  of  particle  diameter  is  shown  in  the 
lower  plot  of  Fig.  8.  The  specific  particle  sizes,  numbers  of  particles,  surface  areas,  and  related  calculations 
for  the  seven  samples  are  given  in  Appendix  A. 

Errors  in  the  estimate  of  the  surface  area  of  the  boron  samples  are  due  to  several  assumptions.  The 
photographs  are  assumed  to  be  representative  of  the  sample  at  large  and  all  particle  sizes  are  assumed  to  be 
equally  visible  in  the  photographs.  The  latter  is  probably  not  the  case  since  larger  particles  were  probably 
preferentially  counted  because  they  obscure  smaller  particles.  I  also  assumed  the  particles  to  be  perfect 
spheres,  which  led  to  a  lower  estimated  surface  area  than  the  actual  surface  area.  This  assumption  is  a  poor 
one  if  the  sample  contains  grossly  irregular  particles.  For  example,  if  the  shape  is  cubic  instead  of  spherical, 
the  ratio  of  the  true  surface  area  to  the  calculated  surface  area  will  be  1.24.  This  assumes  the  volumes  of  the 
square  and  sphere  are  the  same.  The  ratio  can  be  as  high  as  2.3  if  the  diameter  of  the  sphere  is  assumed  to  be 
the  side  length  of  the  cube. 

Comparison  of  Estimated  and  Measured  Surface  Areas 

Table  1  contains  the  measured  and  estimated  surface  areas  for  the  samples  in  this  study.  The  boron 
samples  are  listed  in  order  of  decreasing  BET  measured  surface  area  in  column  5.  Two  different  diameters 
are  listed  in  Table  2.  The  surface  diameter,  d^,  corresponds  to  the  total  surface  area  divided  by  the  total 
particle  number.  The  Sauter  diameter,  d32,  corresponds  to  the  ratio  of  the  sum  of  the  surface  area  to  the  sum 
of  the  mass  of  the  particles.23  Therefore,  one  is  normalized  by  the  particle  number  and  the  other  by  the  particle 
mass  or  volume.  Appendix  A  presents  the  details  of  how  the  diameters  are  obtained. 

The  calculated  surface  area  is  expected  to  be  somewhat  lower  than  the  measured  surface  area  due  to 
the  spherical  particle  assumption  discussed  previously.  In  the  case  of  the  first  two  amorphous  samples,  which 
resemble  spheres,  the  differences  are  expected  to  be  low.  The  calculated  values  are  only  1.3  and  1.2  times  that 
measured  for  samples  1  and  2,  respectively.  Sample  3  has  a  higher  value  of  3.7  and  this  is  partially  attributed 
to  some  nonspherical  characteristics  in  the  SEM  but  primarily  due  to  a  larger  distribution  of  particle  sizes  and 
errors  associated  with  counting  the  particles.  The  SEM  reveals  what  appears  to  be  clustering,  making  it  more 
difficult  to  count  the  particles  or  determine  their  actual  size.  The  crystalline  samples  exhibit  wide  distributions 
of  particle  size  and  vary  considerably  from  spherical.  A  correction  from  spherical  to  cubical  would  yield 
calculated  surfaces  about  1.7  times  higher.  In  the  case  of  sample  6,  however,  it  would  appear  that  the 
difference  is  due  to  a  miscounted  distribution  since  the  sieved  sample,  which  is  more  homogeneous,  is  very 
close  to  the  measured  surface  area.  On  the  other  hand,  the  act  of  sieving  may  select  rounder  particles  and  so 
we  are  not  able  to  clearly  attribute  the  underestimation  of  the  calculated  surface  areas  to  one  particular  cause. 

Sample  4  has  the  worst  match  between  measured  and  calculated  surface  area.  It  may  have  a  more 
rugged  surface,  however  a  slight  miscounting  of  very  large  particles  can  strongly  affect  the  calculated  surface 
area.  Samples  4  and  5  have  similar  distributions  in  Table  2  and  in  Figs.  4  and  5,  however  the  measured  BET 
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surface  areas  differ  significantly,  4.2  and  1.6  m2/g,  respectively.  This  demonstrates  the  importance  of  the 
surface  area  measurement. 

Sample  7's  distribution  was  calculated  first  by  assuming  a  single  value  for  the  size  of  the  particles 
and,  second,  by  counting  the  distribution  of  the  particles.  Both  are  relatively  accurate.  The  distribution  is 
slightly  underestimated.  The  third  contribution  to  the  errors  in  these  calculations  may  be  a  bias  from  not  fully 
counting  smaller  sized  particles.  This  error  would  be  larger  for  larger  distributions  of  particle  size. 

CONCLUSIONS 

Several  conclusions  can  be  drawn  from  comparisons  between  calculated  and  measured  surface  area. 
When  there  is  a  wide  spread  of  particle  sizes,  even  a  counting  method  based  on  SEM  photos  can  have 
moderate  errors.  An  estimate  of  the  surface  area  from  the  radius  of  the  particle  yields  a  minimum  surface  area 
that  can  be  4  to  5  times  lower  than  obtained  from  the  BET  measurement.  The  BET  measured  surface  area 
provides  an  indisputable  comparison  of  the  surface  area  available  for  each  of  the  samples  and  it  accentuates 
the  difficulty  in  getting  accurate  particle  size  distributions.  This  experiment  also  demonstrates  the  difficulty 
in  obtaining  well-characterized,  narrowly  distributed  particles  for  combustion  experiments.  For  example,  even 
after  careful  sieving,  smaller  particles  are  shown  to  adhere  to  the  surface  of  samples  6  and  7,  as  previously 
reported  by  Foelsche  et  al.8 

All  the  crystalline  samples  examined  here  have  larger  particle  sizes  and  lower  surface  areas  than  the 
amorphous  samples.  The  minimum  difference  is  a  factor  of  3  between  samples  3  and  4.  If  combustion 
experiments  determine  that  amorphous  boron  leads  to  more  efficient  combustion,  these  studies  must  be 
followed  by  experiments  comparing  amorphous  and  crystalline  material  in  which  differences  in  particle 
surface  areas  are  eliminated.  If  combustion  experiments  determine  that  crystalline  boron  leads  to  more 
efficient  combustion,  going  to  smaller  crystalline  particles  to  increase  the  surface  area  may  still  improve 
performance.  However,  if  surface  area  and  not  phase  or  purity  is  the  most  important  performance  factor,  the 
amorphous  material  studied  here  should  perform  better  based  on  its  significantly  higher  surface  area. 

Amorphous  boron  is  generally  less  pure  than  crystalline  boron.  However,  sample  1  is  reported  to  be 
99.5%  pure,  which  is  as  pure  as  the  crystalline  samples  are  reported  to  be.  It  may  be  useful  to  have  elemental 
analysis  done  on  the  samples  to  verify  the  reported  purities,  however,  I  did  not  have  enough  material  for  these 
tests.  The  impurities  in  sample  3  are  reported  to  be  primarily  magnesium.13  As  described  by  Neslain,24  the 
particular  impurities  depend  on  the  synthesis  process  used.  In  general,  these  impurities  include  C,  O,  Mg,  Al, 
Si,  Ca,  and  Fe.21  Goodfellow  reports  impurities  of  Fe,  C,  O,  N,  and  H  in  their  99.6%  pure  crystalline  powder 
of  <  1 50  pm  diameter. 

Using  the  model  of  Yeh  and  Kuo,13  the  combustion  time  in  the  first  stage  should  increase 
proportionally  to  the  particle  diameter.  In  the  second  stage  the  combustion  time  should  increase  either 
proportionally  (kinetically  controlled  -  small  particles  at  low  pressure)  or  increase  as  the  square  of  the  particle 
diameter  (diffusion  controlled  -  larger  particles  at  high  pressure).  For  accurate  fitting  of  the  rate  data,  the 
surface  areas  must  be  determined  accurately.  This  report  demonstrates  that  a  direct  surface  area  measurement 
provides  an  accurate  way  of  comparing  different  boron  samples. 


Table  2 


The  0.08  |um  diameter  is  near  the  lower  limit  of  detection. 
Based  on  a  single-valued  diameter  of  25  jam. 
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Foelsche  et  al.8  report  no  difference  in  ignition  times  between  1  pm  and  24  pm  crystalline  boron 
particles  under  high  temperature  (2450  K)  and  high  pressure  (8  to  65  atm)  conditions.  Samples  4  and  7, 
respectively,  were  used  in  their  study.  The  surface  areas  of  these  samples  differ  by  a  factor  of  35.  Others  have 
also  reported  an  insensitivity  to  particle  size.3,5,7 

Rice  and  Russell25  have  observed  more  complete  combustion  products  from  crystalline  boron  (sample 
4)  than  amorphous  boron  (sample  1)  under  high-pressure  (30  Kbar),  laser-initiated  experiments.  The  reactions 
of  boron  with  PETN  and  PETN-NF2  were  compared.  Several  factors  may  play  a  role  in  these  experiments. 
First,  it  is  difficult  if  not  impossible  to  quantify  the  amount  of  boron.  Second,  phase  changes  may  occur  at 
the  higher  pressures.  The  crystalline  sample  clearly  underwent  more  complete  combustion  under  high- 
pressure,  low-temperature  quenching  conditions,  although  the  ciystalline  sample  had  l/9th  the  surface  area 
of  the  amorphous  sample.  These  factors  must  be  kept  constant  or  controlled  to  draw  further  conclusions  from 
this  anecdotal  study. 

Data  and  modeling  from  Yeh  and  Kuo13  taken  with  data  from  Macek3  indicate  a  decrease  in 
combustion  bum  time  with  smaller  particle  size  under  low  pressure  conditions.  It  may  be  that  higher  pressures 
are  the  overriding  factor  in  the  combustion  times  of  the  Foelsche  et  al.8  and  the  Rice  and  Russell25  experiments. 
In  any  case,  better  information  regarding  the  actual  surface  areas  of  the  samples  reported  here  will  hopefully 
clarify  some  of  the  remaining  boron  particle  combustion  issues  surrounding  purity  and  the  effect  of  particle 
size. 
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APPENDIX  A: 

The  following  are  tables  of  the  size  distributions  taken  from  the  SEM  photographs  and  several  other 
calculations  based  on  the  diameter  of  the  particles.  The  two  rows  at  the  bottom  of  the  list  are  the  statistics  for 
the  columns  and  contain  the  average  value  and  the  sum  of  each  column,  respectively.  Column  A  is  the  particle 
diameter.  Column  B  is  the  number  of  particles  in  a  given  area.  Column  C  is  the  interval  size  for  the  particle 
diameter.  Column  D  is  the  surface  area  (SA)  for  a  particle  of  the  corresponding  diameter  in  Column  A. 
Column  E  is  the  surface  area  of  all  particle,  i.e.,  column  D  *  column  B.  Column  F  is  the  mass  of  the  particle 
with  the  diameter  given  in  column  A.  Column  G  is  the  surface  area  per  gram  of  material.  Column  H  is  the 
surface  area  per  particle  times  the  number  of  particles  with  a  given  diameter  normalized  by  the  particle  size 
interval  (column  D*  column  B*  column  C).  Column  I  is  the  normalized  mass  in  grams  obtained  by 
multiplying  the  particle  mass  by  the  number  of  particles  time  the  particle  size  interval  (column  F*  column  B* 
column  C).  The  sum  of  column  J  is  the  numerator  of  the  square  of  the  surface  area  diameter,  D20  , 
corresponding  to  the  particle  diameter  in  column  A.  This  is  defined  in  equation  2: 

(D20  )2  0  =  J  D2  F(D)  dD/  J  D°  F(D)  dD  (3) 

where  F(D)  is  the  particle  size  distribution  function  and  D  is  the  particle  diameter  and  the  integrals  are  taken 
from  zero  to  infinity.  Column  J  is  obtained  by  taking  column  A2  *  column  B  *  column  C.  The  sum  of  the 
denominator  of  this  equation  is  given  in  column  L  and  is  obtained  by  taking  column  B*  column  C.  A  more 
general  equation  is  obtained  by  substituting  j  for  2  and  k  for  0.  The  sum  of  column  K  is  the  numerator  of 
the  Sauter  diameter  equation  given  in  equation  3.  See  reference  23  for  further  discussion. 

(D32  )3  2  =  J  D3  F(D)  dD/  /  D2  F(D)  dD  (4) 

Column  K  is  obtained  by  taking  column  A3  *  column  B  *  column  C.  The  denominator  of  this  equation  is 
the  same  as  the  sum  of  column  J,  i.e.,  the  numerator  of  the  D20  value. 
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